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Spin-orbit coupled superconductors are potentially interesting candidates for realizing topological 
and potentially non-Abelian states with Majorana Fermions. We argue that time-reversal broken 
spin-orbit coupled superconductors generically can be characterized as having sub-gap states that 
are bound to localized non-magnetic impurities. Such bound states, which are referred to as Shiba 
states, can be detected as sharp resonances in the tunneling spectrum of the spin-orbit coupled 
superconductors. The Shiba state resonance can be tuned using a gate-voltage or a magnetic 
field from being at the edge of the gap at zero magnetic fields to crossing zero energy when the 
Zeeman splitting is tuned into the topological superconducting regime. The zero-crossing signifies a 
Fermion parity changing first order quantum phase transition, which is characterized by a Pfaffian 
topological invariant. These zero-crossings of the impurity level can be used to locally characterize 
the topological superconducting state from tunneling experiments. 



Introduction: Majorana Fermions (MF) have been 
the subject of intense recent study, both due to their 
fundamental interest as a new type of particle with non- 
Abelian statistics and their potential application in topo- 
logical quantum computation (TQC)[l|-|5|. Topologi- 
cal superconductors Q are promising candidates for the 
practical solid state realization of MFs [7-14[. A sim- 
ple topological superconducting (TS) system supporting 
MFs, which has attracted considerable experimental at- 
tention p], consists of a spin-orbit coupled semiconduc- 
tor in a magnetic field placed in contact with an ordi- 
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nary superconductor 10l4l4|. It has been shown that the 
semiconductor proposal in one-dimension, i.e. a semicon- 
ducting nanowire, can be driven into a TS phase through 
the appropriate tuning of the semiconductor chemical po- 
tential or equivalently, the carrier density 13|, |14| . Such 



nanowires provide a realization of MFs at its ends in the 



same class as p-wave superconductors [15|, ll6| . The TS 
state is predicted to be realized in the nanowire when 
ever the chemical potential /i with respect to the bot- 
tom of one of the electron sub-bands of the nanowire 
satisfies the constraint < a/A^ + VJ, where Vz is 
the magnetic-field induced Zeeman potential in the wire 
and A is the superconducting pairing potential induced 
in the wire from contact with the s-wave superconduc- 
tor [17| . The s-wave proximity effect on an InAs quantum 
wire, which also has a sizable SO coupling, may have al- 
ready been realized in experiments llSl]. Therefore, the 
TS phase in a semiconductor quantum wire may be one 
of the most experimentally promising approaches to re- 
alizing MFs with non-Abelian statistics. In fact, recent 
experimental results 19] on the semiconducting wire sys- 
tem suggesting the existence of MFs in this system has 
created a great deal of interest in the physics communtiy 

While several schemes such as tunneling and the frac- 
tional Josephson effect have been proposed to detect the 
presence of MFs [13, E, [3, [11, [lH, local probes that 




FIG. 1: (Color) (a) Shiba states in a nanowire (NW), which 
is in a magnetic field B and in proximity to a superconductor 
(SC) can be detected as a conductance peak from an STM 
tip or other tunnel probe near a gate-induced tunable po- 
tential G2. (b) The gate potential produced by G2 can be 
replaced by a local imupurity. (c) Shiba bound state energies 
as a function of inverse impurity strength or impurity trans- 
parency for various values of the Zeeman splitting Vz for 
a topological superconducting nanowire (i.e. Vz = 0.6 meV 
(red curve), Vz ~ 0.8 meV(blue curve)) shows a characteristic 
zero-energy crossings, while the wire in the non-topological 
superconducting phase (i.e. Vz = 0.2 meV (black curve)) 
shows weakly bound Shiba states. 



characterize the TS state of the wire are still missing. 
The need for such local characterization becomes specif- 
ically urgent because of the presence of disorder in real- 
istic semiconducting wires. The absence of time-reversal 
symmetry leaves the proximity-induced superconducting 
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FIG. 2: (Color online) Impurity bound state energy and bulk 
gap as a function of Zeeman splitting for various values of im- 
purity strength g. The values of Vz greater than the critical 
value Vz ~ 0.7 meV, at which the bulk gap vanishes, support 
a topological superconducting phase in the nanowire. The 
Shiba bound state energies are seen to be deeply bound with 
energies signihcantly below the bulk gap only in the topolog- 
ical superconducting phase. Such Shiba bound states might 
contribute to near-zero bias conductance that is robust over 
a range of Zeeman potentials seen in recent experiments [l^ . 
Here we have chosen ^ = 0.2 meV, A = 0.5 meV and the spin- 
orbit coupling a = 0.3 corresponding to a spin-orbit energy 
Eso « SOAteV 



gap of the nanowire susceptible to disorder [23l - l27l | , which 
is predicted to lead to subgap states in the bulk of the 
wire 122 



23[. Such disorder induced gapless states can 



broaden the MFs into Majorana resonances and eliminate 
the fractional Josephson effect [2^, [2^. Therefore, while 
probes of topological superconductivity such as the frac- 
tional Josephson effect, provide a true characterization 
of the topological degeneracy associated with MFs, they 
cannot distinguish if some part of the disordered wire is 
locally in an essentially TS phase. Scanning tunneling 
microscopy (STM) of superconductors has provided one 
such route to experimentally characterizing, in a local 
way, the superconducting state of superconductors such 
as the high-Tc cuprate superconductors 130j. In fact, the 
STM spectrum of impurities allows the characterization 
of not only the superconductor [sij, but also the effect 
of different impurities on the superconducting state [s^ . 
Magnetic impurities in spin-singlet superconductors lead 
to Shiba states, which appear as sharp features in tunnel- 
ing spectra [33| . Tuning the strength of the impurity has 
been predicted, in principle, to lead to a local quantum 
phase transition(QPT) [3J, |35|. Therefore, it is natural 
to expect that studying the tunneling spectra around im- 
purities in wires might allow one to understand the TS 
character of each part of the wire in a local way. 

In this paper, we show that the three ingredients of 
spin-orbit coupling, Zeeman splitting and s-wave su- 
perconductivity required to realize the TS phase in 
nanowires are also the ingredients that realize sub-gap 



Shiba states bound to non-magnetic impurities. Such 
Shiba states in semiconductor nanowires in proximity to 
a superconductor can be probed by a STM or tunnel 
probe arrangement shown in Fig.[IJa,b). For static impu- 
rities, we find that the energy associated with the Shiba- 
state can be tuned to cross zero as seen in Fig.[TJc), pro- 
viding a realization of the previously predicted Fermion- 
parity changing transition [35| in a scenario where it 
might be easier to tune the impurity strength. Such a 
zero-crossing will be shown to be associated with a lo- 
cal version of the Pfaffian topological invariant that may 
be used to characterize the TS phase. By a numerical 
study of models for the semiconducting nanowires, we 
find (see Fig. [5]) that such strongly bound Shiba states 
only occur in the TS phase of the nanowires. The energy 
of these sub-gap states in the TS phase are found to be 
independent of the Zeeman energy and provide a possi- 
ble alternative explanation to the recently observed zero- 
bias conduction peaks [2c|. Such sub-gap states may be 
viewed as a way to destroy the TS gap, which is compli- 
mentary to the Born approximation approach applicable 
to short-scale disorder such as surface roughness of the 
nanowire. 

Local impurity bound states in superconductors: 

To understand the spectrum associated with impurities, 
we consider the Green function Grcr-r'ry'{t",'i''',E) of the 
superconducting nanowire in the Nambu spinor nota- 
tion. Here r,r' represent spatial coordinates on the 
nanowire, cr, a' refer to the spin indices of the electron 
and T, r' = 0, 1 represent the particle- hole index, which is 
needed to describe both the normal and anomalous parts 
of the Green function of the superconducting nanowire. 
The Green function G matrix for an impurity in a super- 
conductor can be calculated using the Dyson equation 

G{rr';E)=Go{rr';E) 

+ J driG^°\rri-E)V{ri)T,G{r^r'-E) (1) 

where V{r) oc gd{r) is the localized impurity potential 
and G^"'' , ,(r:r':E) is the Green function of a clean 
superconducting nanowire with a Bogoliubov de-Gennes 
(BdG) Hamiltonian 



(2) 



In the case of a single- channel wire, ^ represents the 
chemical potential, a is the strength of Rashba spin-orbit 
coupling, Vz the magnetic field induced Zeeman potential 
and A represents the proximity induced superconducting 
pairing potential [ifl, The impurity strength g is 
inversely related to the transparency Z of the impurity 
potential. The matrices Uz, t^, and are Pauli ma- 
trices associated with the indices tr and r respectively. 
The BdG Hamiltonian Hncan also be used to represent 
multi-band wires 22, [s^,!!^ if A*, Vz, A and g are taken 
to be matrices indexed by the channel index. For local 
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impurity potentials, the energy levels of bound states are 
determined from the r = r' = part of the Dyson equa- 
tion Eq. [TJ which is written as 

G(00; ^) = (1 - gG^°^ (00; E)t,)-^G'^"^ (00; E). (3) 

The Shiba bound state appears as a pole in the Green 
function G, which corresponds to a zero of the matrix 



Det[g-^ -G'^^^ {00; E)t,]^0. 



(4) 



the 



Consistent with Anderson's theorem [24 
above equation describing sub-gap states bound to non- 
magnetic impurities, is found to have no solutions in the 
absence of a Zeeman potential i.e. for Vz = 0. The ab- 
sence of Shiba states continues to hold, even in the pres- 
ence of a spatially uniform Zeeman potential Vz 0, 
if the spin-orbit coupling a vanishes, since the Zeeman 
splitting does not affect the wave-functions of the BdG 
Hamiltonian. Therefore, non-magnetic impurities can 
lead to localized Shiba bound states in wires only in 
the presence of a combination of Zeeman splitting and 
spin-orbit coupling, which are precisely the conditions to 
realize a TS phase. 

In contrast, for finite spin-orbit coupling a, Zeeman 
potential Vz and superconducting pairing potential A, 
numerical solutions of Eq. |4] plotted in Fig. [1] show sub- 
gap Shiba states bound to even non-magnetic impurities. 
While bound Shiba states are found to exist in the entire 
range of Zeeman potential Vz > 0, one observes that 
the Shiba states occur deep inside the bulk gap of the 
nanowire, which is also plotted in Fig. [21 only on the TS 
phase of the nanowire, i.e. Vz > Vz,c = 0.5 meV. The 
TS state of the nanowire can be seen identified in Fig. [2l 
by the vanishing of the bulk gap in the wire at Vz = Vz.c- 
Thus, the existence of deeply bound Shiba states may be 
considered as suggestive of the nanowire being in the TS 
phase. 

Shiba states in the TS phase: While the detec- 
tion of a deeply bound Shiba state might suggest the 
nanowire being in a TS phase, it cannot be taken as a 
confirmatory test. Here we show that the evolution of the 
Shiba bound state energy as a function of the strength of 
the impurity g (or magnetic field or density; see discus- 
sion below) may be used as a precise characterization of 
the TS phase. This is seen from the numerical results for 
the Shiba bound state energy in the TS and NTS phases 
of the nanowire, which are plotted in Fig. [T] The Shiba 
bound state energy in the TS phase in Fig. [T] shows a 
characteristic pair of crossings of i? = at both zero and 
non-zero values of the impurity transparency g~^. 

Parity-changing QPT: The crossing of the Shiba 
bound state at E' = at vanishing transparency g^^ ^ 
is characteristic of the TS phase. This follows from the 
fact that an impurity with vanishing tranparency g^^ ^ 
splits the nanowire into two nanowires with a pair of zero- 
energy MFs at the impurity. The pair of zero-energy 



MF split linearly by tunneling across the as the impu- 
rity is tuned away from vanishing transparency. The 
zero-crossing of the energy of a Shiba-state is accompa- 
nied by a change of the Fermion-parity of the ground 
state. This is because the positive and negative energy 
eigenvalues, ei(g^^) > and —ei{g^^) < of the BdG 
Hamiltonian refer to the different occupancies of the as- 
sociated Fermionic state 1. The ground state of the sys- 
tem at ^ corresponds to the positive energy eigenvalue 
being unoccupied. When the impurity strength g~^ is 
tuned from g^^ to g^^, so that the positive and neg- 
ative energy states cross zero and are interchanged i.e. 
^1(92^) < and —ei{g2^) > Oi the system reaches a 
state where the positive energy state 1 is occupied. Such 
a state with a positive energy eigenvalue occupied is an 
excited state of the Hamiltonian with g^^. However, this 
excited state has the same Fermion parity as the ground 
state at g^^ , i.e. before the energy crossing. Since the 
excited state obtained by tuning the impurity to is 
related to the ground state at by adding a Fermion, 
the ground state at g2^ must have a different Fermion 
parity from the ground state at gi^ ■ Therefore, the zero- 
energy crossing of the Shiba state is associated with a 
QPT where the ground state of the system changes its 
Fermion parity. This is analogous to the QPT proposed 
for tunable mag netic impurities in conventional super- 



conductors [35|. The spin-orbit coupled nanowire pro- 



vides a realization of this interesting transition using a 
non-magnetic impurity, which can be tuned by a local 
gate voltage. 

The Fermion parity change associated with zero-energy 
crossings of Shiba states at vanishing impurity trans- 
parency g^^ ^ in the TS phase requires the existence 
of an odd number of zero-energy crossings at finite im- 
purity transparency g~^ 7^ 0. This follows from the fact 
that the total number of zero-energy crossings going from 
g~^ = —00 to g~^ = 00 must be even, since both these 
points are associated with the ground state Hamiltonian 
of the nanowire with no impurity i.e. at 17 ~ 0. Thus, the 
TS phase is characterized by the Shiba bound state en- 
ergy crossing zero an odd number of times as the limit of 
infinite impurity strength — ^ cx) is approached from at 
least one of the sides of either strong repulsive impurities 
or strong attractive impurities. 

Local Pfaffian topological invariants associated 
with impurities: To strengthen the argument that the 
zero-energy crossing is a topological QPT, , we com- 
pute the fermion parity using the Pfaffian topological in- 
variant [3], which can be written in terms of the BdG 
Hamiltonian [3§| . The Pfaffian topological invariant for a 
particle-hole symmetric BdG Hamiltonian HbcIG can be 
written in terms of the particle-hole matrix A = (JyTy 
as QiHsdc) = sgn(P/[7?BdGA]). The particle-hole 
matrix A is defined so that the particle-hole symme- 
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try of Hbclg can be written as HbcLG^ 
which is equivalent to the condition that HbcIG^ is anti- 



SdG' 
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FIG. 3: (Color online) The evolution of the Shiba bound state 
energy with the local chemical potential show an odd number 
of crossings in the topological superconducting phase of the 
nanowire. The chemical potential in the neighborhood of an 
impurity in a nanowire can be varied, in principle, by a local 
gate voltage. The zero-energy crossings occur at different val- 
ues of the local chemical potential for different values of the 
impurity transparency g~^. 



symmetric, which in turn allows the definition of the Pfaf- 
fian PflHEdoM- Defining the QPT of the impurity prob- 
lem requires a local definition of the PfafRan invariant. 
To obtain such a local definition we note that zero-energy 
crossings at an impurity also imply zero-energy crossings 
of the Green function U = G'^"\0, 0; E = 0) - g-'^r, (see 
Eq. |4]) . Since we are considering an impurity in a gapped 
superconductor G^^"^ and the relevant matrix for deter- 
mining zero- modes [/ = G(o)(0,0;£; = 0) — g is both 
Hermitean and particle-hole symmetric. Therefore the 
local Pfaffian invariant is defined as 

Q(g-i) = sgn (P/[(G(0) (0, 0; E ^ 0) - ff" V,)A]) . (5) 

The above Pfaffian invariant can only change sign when 
the determinant of J7 = G'^°'> {0,0; E = 0) - c^'V^ van- 
ishes, which corresponds to a solution of Eq. \^a.t E — 
and therefore a zero-crossing of Shiba bound state en- 
ergies. The Fermion parity changing transition of the 
ground state as a function of impurity strength g is char- 
acterized by the topological invariant Q{g~^) changing 
sign. Therefore, the local Fermion parity at the impu- 
rity can be calculated directly from computing the lo- 
cal Fermion-parity Q{g~^) without computing the Shiba 
bound states. 

Even though the strength of an electrostatically in- 
duced impurity in a wire can in principle be controlled 
by a gate voltage in a simpler way than the tunable rnag- 
netic impurity required for s-wave superconductors [3a |. 
tunable non-magnetic impurities also might be difficult 
to obtain. An alternative approach to locally character- 
izing a TS phase is to study the evolution of the Shiba 
bound state energy as one tunes the chemical potential 
or the magnetic field towards the gap-closing topological 



phase transition. It is reasonable to expect a sufficiently 
strong impurity behaves has an asymptotically stronger 
effect as one approaches the phase transition. The evo- 
lution of the Shiba states across the phase transition for 
different local impurity strengths, which are plotted in 
Fig. [21 dotted green) and [31 shows that this is indeed 
the case. Therefore, the TS phase of the nanowire can 
also be characterized by the presence of an odd number 
of zero-energy crossings of the Shiba bound state energy 
as the applied chemical potential is locally tuned across 
a QPT. 

Conclusion: Spin-orbit coupled coupled nanowires 
together with Zeeman splitting and s-wave superconduc- 
tivity, which can be used to realize one-dimensional TS 
nanowires, are also precisely the ingredients needed to 
realize sub-gap Shiba states bound to non-magnetic im- 
purities. Similar to magnetic impurities s-wave super- 
conductors [3l|, Shiba bound states can also be used 
to characterize the TS phase of the nanowire. Specif- 
ically, only the TS phase of the nanowire is found to 
support deeply bound Shiba states, whose energies are 
only weakly Zeeman field dependent. Such low-energy 
Shiba states found only in the TS phase provides a pos- 
sible alternative explanation to the recent observations 
of zero-bias conductance peak, which was interpreted as 
a signature for the existance of MFs [3]. In the TS 
phase of the nanowire that the Shiba bound state energy 
is found to cross zero-energy an odd number of times 
as one tunes the strength of the impurity. Such zero- 
crossings of the Shiba bound state energy, besides pro- 
viding a local characterization of the TS phase of the 
nanowire, are also intrinsically interesting since they are 
associated with a Fermion parity changing QPT of the 
nanowire ground state (35| . For systems where the im- 
purity strength is difficult to tune, we find that similar 
zero-energy crossings of the Shiba bound state energies 
for impurities in nanowires in the TS phase can also be 
obtained by gate tuning the chemical potential. There- 
fore, the TS phase of spin-orbit coupled semiconducting 
nanowires can be locally characterized by studying the 
evolution of Shiba bound state energies, which appear as 
sub-gap resonances in the tunneling spectrum into the 
wire. Since most of our results follow from the scattering 
equation (Eq. [3]) from a local impurity, we expect our 
conclusions to apply to multi-channel nanowires as well 
as impurities in two-dimensional TS systems. Such local 
characterizations of the TS phase of the nanowire are par- 
ticularly useful, since only parts of wires can be expected 
to enter the TS phase in disordered semiconducting wire. 
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